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It is shown that finite size effects in the free energy of a rough interface of the 3D Ising and three-state Potts 
models are well described by the capillary wave model at two-loop order. The agreement between theoretical 
predictions and Monte Carlo simulations strongly supports the idea of the universality of this description of 
order-order interfaces in 3D statistical systems above the roughening temperature. 



1. Introduction 

Many 3D physical systems show the presence of 
interfaces separating coexisting phases at thermal 
equilibrium. It turns out that these interfaces are 
often dominated by long wavelength fluctuations 
(i.e. they behave as fluid interfaces). These soft 
modes play an essential role in the description of 
finite size effects (FSE) in the fluid interface's free 
energy (see e.g. [1]). 

3D spin systems offer a simple context where 
these effects appear and can be studied, by using 
numerical simulations to check theoretical predic- 
tions. It is in fact well known that at low tem- 
perature, on finite volumes, they show domain 
walls separating coexisting phases, which behave 
as fluid interfaces between the critical and the 
roughening temperature. Besides, the models we 
study are interesting because they are related by 
duality to 3D Z^ gauge theories, and to high tem- 
perature 4D gauge theories through dimensional 
reduction. 

While below the roughening temperature the 
interfaces are almost rigid and a theoretical mi- 
croscopical approach can be taken (see e.g. [2] 
and references therein), above it one is forced to 
assume an effective model describing the collec- 
tive degrees of freedom of the rough interfaces. 

We follow the capillary wave model (CWM) [3], 
in its simplest formulation, assuming an effective 



hamiltonian proportional to the area A[x] of the 
surface 

where the single-valued function x(r,t) describes 
the displacement from the equilibrium position of 
the interface and a is the reduced (order-order) 
interface tension. To compare the predictions of 
the CWM with numerical results obtained from 
Monte Carlo (MC) simulations, one can choose 
3D lattices of R x T x L sites, with L ^> R,T and 
periodic boundary conditions in each direction. 
This particular choice of the lattice shape allows 
one to consider only interfaces orthogonal to the 
elongated direction L, the probability of having 
interfaces orthogonal to the other directions be- 
ing negligible. The 2D field x(r,t) is therefore de- 
fined on the rectangle (r,t) £ [0,R] x [0,T] with 
opposite edges identified, i.e. on a torus. 

2. The two-loop CWM approximation 

Rather strong FSE, depending on the shape of 
the lattice, already appear at the one-loop (gaus- 
sian) approximation to the CWM [4]. The par- 
tition function Z cw = J [Dx] e ~ aA ^ can be in 
fact expressed as an expansion in powers of the 
adimensional parameter cr" 1 , proportional to the 
minimal area of the surface (cr a = aRT), 
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Zcw{z, a a) = Se a " ■ Z(u)(z) ■ Z(2i){z,o- a ) ... (2) 
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where S is an unknown constant and z = R/T. 

The one-loop contribution ^(i/), obtained re- 
taining only the quadratic term in the expan- 
sion of A[x], depends only on z, namely on the 
asymmetry of the transverse sizes of the elon- 
gated lattice. is nothing but the (exact) 
partition function of a 2D conformal invariant free 
boson on a torus of modular parameter r = iz: 
Z (U)( Z ) = 7? \V ( iz ) h 0')l~ 2 > where 77 is the 
Dedekind function tj(t) = q 1 / 24 Yl^ =1 (1 — q n ) , 
q = e 27 ™ r (gee e.g. [4] and references therein). 
The one-loop CWM prediction has been tested, 
on asymmetric lattices (R ^ T), in the scaling 
region of the 3D Ising model and a quite remark- 
able agreement with the MC data has been found 

Taking into account higher order corrections to 
the gaussian model a more stringent test on the 
CWM can be provided: in fact, while many differ- 
ent effective hamiltonians reduce to the gaussian 
form at one-loop level [5] , they differ in the form 
of two- and higher-loop corrections. 

The two-loop term Z^i) of Eq.(2) can be calcu- 
lated perturbatively, expanding (1) at the next- 
to-leading order in cr a , and is given by [6] 

Z i2l) (z,a a ) = l + -^- f(z) + 0{a-J) , (3) 

where 

/(z)^{[^e 2 (7z)] 2 -£ze 2 (iz) + !} (4) 

and e 2 (r) = 1 - 24£~ =1 f^, q = e 2 " T , is the 
second Eisenstein series. 

In contrast to what happens at the gaussian 
level, here the capillary wave contributions do not 
depend only on the asymmetry parameter z but 
also on a a . As a consequence, the two-loop finite 
size behavior of the energy splitting E occurring 
between vacua on finite volumes 

E(R,T) = Z cw (z,a a ) (5) 

has no longer, the well known classical functional 
form [8] E d (R) = 8e~ aR not even for symmetric 
(R = T) lattices, as can be easily seen by setting 
z = 1 into Eq.(5). 



3. Numerical results 

The comparison between formula (5) and the 
values of E extracted from MC simulations pro- 
vides a simple and stringent way to verify the 
CWM predictions. Notice that no new free pa- 
rameters are introduced within this approach: 
Eq.(2) contains the same number of undeter- 
mined parameters, namely a and S, to all orders 
of approximation. 

We report here the main results of the MC sim- 
ulations we have made for the three-state Potts 
[7] and the Ising model [9] : in both cases we have 
followed Ref.[10] to extract the energy splittings 
E. For the Potts model, in order to deal only 
with order-order interfaces, the simulations have 
been performed at [3 = 0.3680 (in the conventions 
in which f3 c ~ 0.36708). The results are plotted 
in Fig.l: the sizes vary in the range 9 < T < 20, 
10 < R < 36, with R > T, and the longest lattice 
size is fixed at L = 120. 
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Figure 1. CWM at two-loop with the MC data 
for the Potts model at [3 = 0.3680. The lines 
represent the best fit of all data to Eq.(5), a and 
S being given in Tab.l. Error bars not reported 
are smaller than the plotted symbols. 

For the Ising model, the simulations have been 
performed at /? = 0.2246 and /? = 0.2258. The 
results are plotted in Fig. 2. It is evident that no 
single valued function of cr a can describe the data 
reported in both figures. In particular, within the 
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Table 1 



Comparison of interface tensions obtained from different theoretical approximations. 





approx. 


a 


S 


xVd.o.f. 


Potts 0.3680 


2-loop 


0.00991(7) 


0.138(2) 


0.73 




1-loop 


0.01005(7) 


0.172(2) 


3.60 




class. 


0.00809(7) 


0.139(2) 


36.3 


Ising 0.2246 


2-loop 


0.00649(10) 


0.108(2) 


0.58 




1-loop 


0.00693(10) 


0.142(3) 


1.65 




class. 


0.00651(10) 


0.131(3) 


11.1 


0.2258 


2-loop 


0.00941(6) 


0.123(2) 


0.44 




1-loop 


0.00958(6) 


0.148(2) 


1.48 




class. 


0.00833(6) 


0.121(2) 


19.2 



classical approximation (Z cw = Se~ a ") of Eq.(5) 
to fit each sample of data, at fixed [3, one always 
finds a bad \ 2 /d.o.f.. A clear enhancement of 
the latter can instead be seen if one uses the one- 
loop approximation (Z cw = Se~ a " Z^{z)), while 
good C.L. are reached with the full two-loop form 
ofEq.(5). 
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Figure 2. The same as Fig.l in the case the 
Ising model at /? = 0.2246 and /? = 0.2258. 

The results of these fits are given in Tab.l. 
Thus, the CWM in the two-loop approximation 
provides an excellent description of order-order 
interfaces in both the Ising and the three-state 
Potts models. This result is a strong indication 
of the universality of this description of interface 
physics in 3D statistical models. 
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